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Abstract 



We improve the calculation of the supersymmetric 0(a s ) QCD corrections to 
the decays of Higgs bosons into quarks and squarks in the Minimal Supersymmetric 
Standard Model. In the on-shell renormalization scheme these corrections can be 
very large, which makes the perturbative expansion unreliable. This is especially 
serious for decays into bottom quarks and squarks for large tan/3. Their corrected 
widths can even become negative. We show that this problem can be solved by a 
careful choice of the tree-level Higgs boson couplings to quarks and squarks, in terms 
of the QCD and SUSY QCD running quark masses, running trilinear couplings A q , 
and on-shell left -right mixing angles of squarks. We also present numerical results 
for the corrected partial decay widths for the large tan (3 case. 



1 Introduction 



The Minimal Supersymmetric Standard Model (MSSM) [p]] is a very promising extension 
of the Standard Model. This model has five physical Higgs scalars, h°, H°, A , and H ± . 
Studying their properties is very important for probing the nature of the breaking of the 
electroweak symmetry and supersymmetry. 

The Higgs scalars have large couplings to quarks and scalar quarks (squarks) of the 
third generation 0. Therefore, in typical cases the decays to top and bottom quarks 
have large branching ratios and are often the main modes, see e.g. H . Decays to top 
and bottom squarks can be also dominant [|], H if they are kinematically allowed and the 
left-right mixing parameters of the squarks are large. Studying these decays is therefore 
very important for the discovery and the detailed study of the Higgs bosons in the MSSM. 

The decays of Higgs bosons into quarks and squarks receive large radiative corrections 
by QCD interactions. In the MSSM, one-loop SUSY QCD corrections to the decay widths 
have been given analytically, both for quark modes |6|, [7|, ||, Q and squark modes jR], |TT 



However, when the lowest order widths are given in terms of the on-shell parameters for 
quarks and squarks, the 0(a s ) correction terms are often comparable to or even larger 
than the lowest order ones, especially for decays into bottom quark and squark. Such 
large corrections make the perturbation calculation of the decay widths quite unreliable. 
When the correction term is negative, the corrected width can even become negative, 
which clearly makes no sense. This "negative width" problem has also been observed in 
other processes of squarks ||12 . 



The couplings of Higgs bosons to quarks and squarks depend on the quark masses m q . 
For the decays into quarks, it is well known |L3|, [14]] that the large part of the gluon loop 
correction can be absorbed by giving the tree-level Higgs-quark couplings in terms of the 
QCD running quark masses m g (Q) SM , where Q is at the scale of the parent Higgs boson 
mass. The convergence of the perturbation series is greatly improved by this replacement. 
However, the gluino loop correction to the decays into bottom quarks can also be very 
large for large tan (3, the ratio of the vacuum expectation values of two Higgs fields. It 
has been pointed out [15| that the main source of this correction is the counterterm for 



m&. Its contribution to the decay widths into bottom quarks was included in ref. |T6[ by 
effective Higgs-bottom couplings with one free parameter. However, no improvement of 
the large gluino loop corrections has been presented along this line. Moreover, the Higgs 
boson couplings to squark are functions of m q , Higgs-squark trilinear couplings A q , and 
squark left-right mixing angles 9q. The corresponding improvement of large corrections 
to squark modes is therefore more complicated and has not been discussed so far. 

In this paper we improve the one-loop SUSY QCD corrected widths of the MSSM 
Higgs boson decays into quarks and squarks. We concentrate on the decays into bottom 
quarks and squarks, where the large QCD corrections are most serious in a phenomeno- 
logical study. The essential point of the improvement is to define appropriate tree-level 
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couplings of the Higgs bosons to quarks and squarks, in terms of the running quark 
masses m q (Q) both in non-SUSY and SUSY QCD, running Higgs-squark trilinear cou- 
plings A g (Q), and on-shell left-right mixing angle 9 q of squarks. Such a treatment of 
the Higgs boson couplings to quarks and squarks will also be useful in studying radiative 
corrections to other processes of Higgs bosons and squarks. 

This paper is organized as follows. In section 2, we review the improvement of the 
large gluon loop corrections to the decays into quarks by using running quark masses. 
In section 3, we discuss the origin of the large gluino loop corrections to the decays into 
bottom quarks in the large tan/3 case, and work out a suitable method of improvement. 
The corresponding improvement of the squark modes is discussed in section 4. Section 5 
shows numerical results of the partial decay widths of Higgs bosons into bottom quarks 
and squarks. Section 6 is devoted to conclusions. Throughout this paper we adopt 



notations and conventions of 10 



2 Gluon loop corrections to Higgs decays to quarks 

We first review the method to improve the calculation of the gluon loop corrections to 
Higgs decays to quarks — > qq (0 = h°, H°, A ) by using QCD running quark masses 
and renormalization group equation, following |13], [TJJ] . A similar discussion holds for the 
decay — > qq'. 

When the Yukawa couplings in the lowest-order decay widths of a Higgs boson are 
given in terms of the pole masses M q of the quarks, the one-loop QCD corrections to 
the inclusive decay widths (including real gluon radiation) have large 0(a s \n(m^/M q )) 
terms from gluon loops. The explicit form of the corrected width is Jl3|, |i~4H , in the limit 

Of M q <C J710, 



1 OS 



1 + 

7T 




The correction can be very large and causes bad convergence of perturbation series. This 
problem is especially crucial for the decays to b quarks. In calculations using dimensional 
regularization with renormalization scale Q = m^, the above correction comes from the 
counterterm 6m q for the Higgs-quark Yukawa coupling h q oc m q . It is well known |L3|, [14| 
that this correction can be absorbed into the tree-level widths by using the (non-SUSY) 
QCD running quark masses m q (Q ~ tti^sm for the tree-level Yukawa couplings. We can 
also resum large higher order corrections by using renormalization group equations for 
quark masses. 

The one-loop relation between M q and the MS (dimensional regularization with mod- 
ified minimal subtraction) running mass m q (Q)sM in the non-SUSY QCD is 

m q (Q) SM = M q + 5m q 9 \Q). (2) 
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The counterterm 5m^(Q) from the gluon loop is at one-loop level 

h », (q)= _-Mm„(.„|;^). ( 3 ) 

One can see that the large logarithmic part of Eq. ([!]) is cancelled when m q (m^) sm is used 
in the tree-level width. 

Next we rewrite Eq. (0) as 

rn q (Q) SM =( m ^ M )rn q (M q ) SM . (4) 
\m q {M q ) SM ) 

In Eq. (H) a large logarithmic correction appears only in the ratio m g (Q)sM/w 9 (M g )sM- 
We can then resum large higher-order corrections of 0(a™ ln n (m,p/M q )) in Eq. (Q) by 
using the renormalization group equations (RGE) for m q (Q)sM- The one-loop running 
a s (Q) is given by 

12tt 

as(g) -(33-2n / )ln(QVA^)' (5) 
The solutions of the two-loop RGEs JI3 are 



^ { J) (33-2n,)lB#- (33-2n^ mf- 1 ' 



7 v "/ 

and for the running quark mass at M q 



m q (M q ) SM = M, 



q 



1 , 4 «j 2 )(M g ) | K ( af){M q 

3 7T q \ 7T 



(7) 



with Ifj ~ 10.9 and lf 6 ~ 12.4 and n/ = 5 or 6 for M 6 < Q < M t or Q > M t , respectively. 
Using the functions 

12 

/23 \23 

c 5 (x) = (— x) (1 + 1.175a;) (M b < Q < M t ) , (8) 

c 6 (x) = C-xX (1 + 1.398a;) (Q > M t ) , (9) 
the solution for the ratio (m q (Q)sM/ m q(M q ) SM ) is 

M t <Q<M« : , (10) 

Ce(af»(M,)/ir) c 5 (a« 2) (M,)/i) 

For the numerical calculations in section 5 we take Q = m^. For the calculation of 
m q(Q)sM we use eq. @) and eqs. (||) - (JTTJ) . For all other calculations eq. (|5|) is used. 
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3 Gluino corrections to quark Yukawa couplings 



3.1 Correction to ra& 

The gluino loop corrections to the decay widths of — > bb J7J and H + — > tb [|| |9| can 
be very large for the cases with large tan/? = v 2 /v\, large m§, and large Their main 
part comes from the counterterm 5m^ by squark-gluino loops which has an enhancing 
coefficient tan (3. In contrast to the case of gluon loops, however, these corrections are not 
caused by a large logarithm. Rather, as we will see now, it originates from the effective 
bbH 2 coupling which is generated by the soft SUSY breaking in the loops. 



Here we follow the discussion in ref . [ [TB, [L7| . For illustration, let us consider the case 



where squarks are much heavier than Higgs bosons. In this case, the interactions between 
Higgs scalars and quarks are properly described by the effective two Higgs doublets theory, 
after integrating out squarks and gluino. The effective interactions of Higgs scalars and 
b which are proportional to mj are written as 

= -h b b R (b L Hl - tjMi) - h b A b b R (b L H°* + t L H 2 ) + (h.c). (12) 

At the tree-level, the bbH 2 coupling is forbidden by SUSY, namely Aft = 0. However, 
this coupling is generated by loop corrections with soft SUSY breaking. Indeed, an explicit 
calculation of the squark-gluino loops for zero external momenta gives the result 

2a 

A b = --^rn sf JLl(mg,MQ,Mf,), (13) 



where 



I(m 1 ,m 2 ,m 3 ) = C (0, 0, 0;m 1 ,m 2 ,m 3 ) 

m\m\ ln(ml/m 2 ) + m\m\ ln(w|/ m|) + m\m\ ln(mg/ m\) 
(rrii ~ m i)( m 2 ~ m i)( m i — m i) 



(14) 



Co is the standard three-point function ]TSj in the convention of [HJ . We see that A& = 
for the exact SUSY case Trig — 0. As long as (jrig, |/^|) are not much larger than Mq ^, 
Aft itself is smaller than unity and therefore does not destroy the validity of perturbative 
calculation. Nevertheless, as shown below, Aft is responsible for the very large gluino 
correction to the — > bb and H + — > tb decays. 

We then evaluate the couplings of Higgs bosons to bb in the mass eigenstate basis. For 
reference, we list the basic relations between the gauge and mass bases 0: 

^ 13 1 = 1 ™" ~Z a a ](".] + ( I' I • US) 



Hi 




cos a 


H° 2 


)-i 


sin a 


m I 


Hi N 


H 
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(17) 



( — cos (3 sin (3 
\ sin (3 cos (3 

G 0,± are the would-be Nambu-Goldstone bosons and irrelevant here. In the large 
limit, cos a — ► sin/5 and sin a — > — cos/3 follows. 
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In the mass basis, Eq. ([12]) becomes (i> 

= _^[cos/3 + A 6 sin/3]&& 
v2 

h 

y= [cos a + At sin a] H°bb 

v2 

+ — -p [sin a — A b cos a] 
v 2 

+ ^ [sin /3 - A b cos $ A°6 75 6 
v 2 

+h b [sm /3-A b cos /3] H'b R t L + (h.c). (18) 
The first term gives the (non-SUSY) QCD running mass m&|sM at scale Q, 

m b (Q)sM = ^=[cos/3 + A 6 sin/3]. (19) 

The difference from the SUSY QCD running mass mf,(Q)MSSM = h b v cos /3/>/2 in the 
DR scheme (dimensional reduction with modified minimal subtraction) is due to the A b 
term. For tan/3 ^> 1, the effect of A b on m b is enhanced by a large vacuum expectation 
value of H 2 . As a result, the relative correction by squark-gluino loops (m b (Q)su — 
^&(Q)mssm)Atc6(<3)mssm has a factor tan/3 and can become very large, even 0(1). This 
property has been pointed out in [DJ, mainly as a weak-scale threshold correction to the 



bottom-tau Yukawa coupling unification. 

In Eq. (|Tj|) the contributions of A b to the Higgs-bottom couplings take forms different 
from those to m b . When the tree-level couplings are given in terms of m b {Q) sm, the 
corrections by A b can be very enhanced for tan/3 3> 1, as pointed out in [fl7l, [3~6f| . For 
example, A°bb coupling is expressed as 

b ih b< 



a = —= [sin (3 — A b cos /3] 
v2 

= = — tan [3 



1 



sin (3 cos (3 1 + A b tan /3 



(20) 



Similarly, the corrections to the couplings H°bb, h°bb, and H + ib can also be very large. 
We see that A b is the main source of the large squark-gluino loop corrections to decay 
widths of (H°, A ) — > bb JF| and H + — > tb || in the on-shell renormalization scheme. Since 
the largeness of the gluino correction comes from the property that the H^b coupling is 



6 



forbidden at tree-level, higher order corrections are not expected to be larger than the 
one-loop corrections. 

Here we have two comments. First, similar to the H 2 bb coupling, the squark-gluino 
loops also generate an effective H\it coupling with a coefficient A t h t . As consequence, 
for large tan/3 large gluino loop corrections of 0(a s A t tan/3) to the decay widths of 
(A°,H°) — > ti appear in the vertex corrections [[?[]. However, their widths and branching 
ratios decrease as tan/5 increases. The contribution of A t to the width of H + — > tb is 
sufficiently small. Therefore, the phenomenological importance of the A t correction is 
smaller than that of A&. Second, the effective coupling A& can be also generated from 
other loop corrections such as the higgsino loops. This effect has also been discussed 
T5], O, EUl PJ, but it is beyond the scope of this paper. 



3.2 Method of improvement 

It is clear from Eqs. (|l8l)-(f2~0|) that the QCD expansion of the Higgs decay widths to 
bottom quarks around the tree-level ones in terms of Mb or m^m^) sm may cause bad 
convergence. As in the case of the gluon loops, we can improve the QCD perturbative 
expansion by changing the definition of the tree-level coupling of bottom quarks and Higgs 
bosons. Ref. |16|] calculated the decay widths into quarks by using effective couplings in 
Eqs. (fl8f)-(]20D and added the gluon loop corrections. However, numerically, A& was not 
calculated from the SUSY parameters. Here we take a different approach, which is more 
suitable for our purpose to improve the SUSY QCD correction. 

The simplest case is the decay A — > bb. The coupling a b , Eq. (|20D , is expressed in 
terms of the SUSY QCD running mass m b (Q) MSSM at Q = as 

a b = Zmfe(g)MSSM taafll - A 6 cot/5]. (21) 

v 

In contrast to Eq. (0), the correction is now very small. This is due to the property that, 
for tan/5 ^> 1, A is almost ImH® while A& is the coupling to H2. We therefore expect 
that the SUSY QCD running mass m^Q = m^MSSM is an appropriate parameter for the 
tree-level A — > bb decay. Since H ± is also almost for large tan/5, 7B;,(7«#±)mssm is 
appropriate for the H + — > tb decay. 

The case of (H°, h°) decays needs a special treatment. The couplings s b H for H°bb and 
s b h for h°bb in terms of rrib(Q)MSSM receive relative corrections A^tana and — A^cota, 
respectively. One of them can be very large. For example, for very large itla, s b h becomes 

4 - - mfe(g)MSSM [l + A b tan/5] = ~ m&(Q)sM . (22) 

v v 

This is the SM decoupling limit for h°. In this case rrib{Q) sm is more appropriate than 
^6(Q)mssm for the tree-level coupling s b h . 
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To find a better way, we make a change of the basis of the Higgs doublets to 

H sm = cos PH{ + sin (3H 2 , H cx = - sin fiH{ + cos fiH 2 , (23) 

where H\ = (—Hi, Hf*). Only H sm has a vacuum expectation value. H sm can be regarded 
as the "true Higgs", while H ex is the "extra scalar doublet". (A , H^) are included in 
H ex . H° and h° have components of both H sm and H cx . 

The couplings of H sm and H ex are obtained from Eq. (0). One can see that the 
H sm bb coupling takes the form h b cp(l + A&tan/3), which is properly parametrized by 
mb(Q)sM- On the other hand, the H ex bb coupling is h b sp(l — A^cot (3), which is properly 
parametrized by m b (Q) MSSM . Using the relation between (H°, h°) and (H sm , H ex ), 

the appropriate choices for the tree-level couplings are 

b I n\ m b(Q)sM . . ( flU n m b(Q)MSSM / -\ 

s)t = — cos(a — p) = h sm a — p) tanp £ , (25) 

s h = sin (a — p) = h cos(a — p) tanp = . (26) 

Eq. ( p2j) is reproduced in the large limit, using cos(a — (3) — > 0. 

In the numerical calculation, we obtain m;,(Q)MssM from m b (Q)sM in Eq. (|J) by using 

m b (Q) M ssM = m b (Q)sM + Smf\Q), (27) 

with the full one-loop contribution of gluino loops and the conversion term between the 
MS and DR schemes 

a 



M q (Bi(M^,ml,ml) + Bi(M' q ,m\,m\ 2 ) + 1) 
+m § sm2^(S (Mj,m?,m? 1 ) -^(M^m^ml))] , (2S 

instead of the zero momentum approximation Eqs. fll3|), (|T4]), and (|20|) . 



The formulae of the effective couplings in this section were obtained in the limit of 
rriA <C rrig. Nevertheless, we apply these formulae in section 5 to the case of rriq < tua- 
This is justified because for the decay modes discussed in section 5 the modification of A b 
and At by large external momentum does not affect the effective Higgs-quark couplings 
significantly. The situation is different for the decays to ti. We therefore leave such decays 
for future works. 
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4 Higgs— squark couplings 



We next consider the SUSY QCD corrections to the Higgs boson decays to squark pairs. 
As in the case of decays to quarks, an appropriate choice of the tree-level couplings of 
Higgs bosons and squarks is essential for improving the QCD perturbation calculation. 



We first review the on-shell renormalization of the squark sector, following ||10|| . For 
given values of /i and tan f3 the masses and mixings of squarks in a generation are fixed by 
five independent parameters, in addition to the masses of the quark partners. We can use, 
as such parameters, the pole masses (m^, m^ 2 , , fn^) and the on-shell left-right mixing 
angles (% 6~ b ) which are independent of the DR scale Q, with one constraint by SU(2) 
gauge symmetry. The on-shell mixing angle 6q is defined by specifying the counterterm 
S9q = #,j(DR) — 9q (on-shell). Various definitions have been proposed in previous works 



on the squark interactions |2T], [XT], [23|]. Here we adopt the definition in [^3J, where 86 q 
is fixed to absorb the anti-hermitian part of the squark wave-function renormalization: 

l Re{nf 2 K) + Hf 2 K- 2 )} 

86g = - (29) 

2 m| - to| 

with nf 2 (p 2 ) the off-diagonal squark self energy in the squark mass basis. Eq. (^9[) is also 
applicable to loop corrections other than in QCD. 

We can then define the on-shell squark parameters {Mq, Mjy, M^, A t , A^) in terms 
of the above on-shell parameters by using tree-level relations. For example, on-shell A q 
is given by 

I 

MqAq = -(m| - m\ 2 ) sin 26 q + M g /i{cot 0, tan/5} , (30) 
where cot (3 (tan/5) is for q = t (6). Also, on-shell M~(g) for the q±^ sector is given by 

M 3(?) = rnl cos 2 6q + m| sin 2 6~ q - m% cos 2/3(1™ - e q sin 2 6 W ) - M\ (31) 

(with I q L the third component of the weak isospin and e q the charge of the quark q). 
These parameters were used as inputs in JTDj , T2\. Note that the value of on-shell M?(t) 



is different from that of M^(b) by finite QCD corrections [0, |12|, 24 



The QCD correction to the decay width of Higgs bosons into squarks consists of four 
parts: vertex correction, squark wave function correction, counterterm for the Higgs- 
squark coupling, and real gluon radiation. Since a squark couples to both Hi and H 2 at 
the tree-level, the vertex corrections are not expected to give a large correction. A large 



correction to decays into b in the on-shell scheme comes mainly from the third part JT0[ , 
the counterterms for the Higgs-squark couplings. These couplings depend on three pa- 
rameters which receive QCD corrections, m q , 6 q , and A q . In the on-shell renormalization, 
their counterterms can be very large and make the perturbation calculation unreliable. 
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As in the decays to quarks, we can improve the perturbation calculation by us- 
ing SUSY QCD running parameters niglm^MSSM and A q {m^) in the tree-level Higgs- 
squark couplings. However, the mixing angles 9g are kept on-shell in order to cancel the 
(9(l/(m| 2 — 77i| )) term in the off-diagonal squark wave function correction, which causes 
a singularity for m, 2 ~ m^. In addition, the counterterm Eq. (|29[) largely cancels the 
0{m q m~g) term in the off-diagonal squark wave function correction, which is large for 
q = t. 

The renormalization of A b requires special attention: For a given set of (rag , mg , 
Mf,), the value of the running A b can be numerically very different from the on-shell 
A b , especially for large tan/3. In Eq. (|30|) there are two sources for a large difference 
5A b . One is the correction due to large 5m b . The other is the correction to the first 
term of Eq. (|30"D, which is the left-right mixing mass of b, M\ R = m b (A b — /itan/?), at 
the tree-level. The QCD correction to this term is of the order of a s M\ R and gives a 
contribution 8 A b = 0(a s /j,taja(3), which is very large for large tan /3. The on-shell A b is 
therefore physically inconvenient, at least for the large tan/3 case. In fact, this is due to 
the property that 9~ b is insensitive to A b . 

Instead of on-shell A b , we can use the on-shell 6~ b as an input parameter. However, 
in order to avoid fine tuning of the parameters so that we do not get a too large value of 
running A b , we rather use A b (Q = m^) as input. 

5 Numerical analysis 

In this section we show the full one-loop SUSY QCD corrected widths of the MSSM 
Higgs boson decays into bottom quarks and squarks (including real gluon emission to 
avoid infrared divergence), namely (h°,H°,A°) — > bb, H + — > tb, (H°,A°) — > b^b*, and 
H + — > t,jb* = 1, 2), with and without the improvement worked out in this paper. 

5.1 Procedure of the numerical calculation 

Our input parameters are all on-shell except A b which is running as explained above, i. e. 
we have M t , M b , Mg(t), M^, M^, A t , A b (Q), fi, tan/3, m^, and nig, with Q at the mass of 
the decaying Higgs boson. (Since real gluinos do not appear in our processes we take 
rrig as a tree-level input, neglecting radiative corrections to m g .) For the kinemat- 
ics (phase space) the on-shell masses are used. For consistency, all arguments in the 
Passarino-Veltman integrals, which appear in the calculation of the one-loop corrections, 
are also taken on-shell. According to the procedure of DR renormalization we set the UV 
divergence parameter A (= - — 7 + log47r) to zero. In the following we will describe in 
detail the procedure for obtaining all necessary on-shell and DR parameters for quarks 
and squarks, namely running m q) pole m 9 - 12 , running (Mq, M^, M^, A t ), and running 
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and on-shell 9$. 



Top— stop sector: The on-shell masses , my 2 , and the mixing angle 9$ are calcu- 
lated by diagonalizing the well-known stop mass matrix in the t^-tR basis. Using Eqs. ([3]) 
and (D we get m t (Q)sM, and using Eq. fl28|) m t (Q) = m t (Q) M ssM = m t {Q)su + Sm[ 9) 



H 

The counterterms Sm~ follow from 

5m\ = Re[n?/ 9) (m|) + + nj®] , (32) 

using Eqs. (25)-(27) of |nj. For 5% we take Eq. (||) with IT? 2 (p 2 ) as given in 0. Note 



that the self-energies 11^ depend on the DR scale Q. Now we can compute the running 
stop masses m~ (Q) = m~ + 5m~ and 9f(Q) = 9%+ 59 f Inserting the running masses and 
mixing angle into the formulas 

M| = m| cos 2 9~ t + m\ sin 2 9~ t - m\ - D L (t) , (33) 

M 2 G = m\ sin 2 0f+ m| cos 2 9 t - m 2 t - D R (t) , (34) 

m t A t = (m| — m| 2 ) sin ^ cos 9t + m t fi cot /? , (35) 

with the Z)-terms 

Dl{q) = m \ cos 2/3 (I q L — e q sin 2 #iy) , D fl (g) = m| cos 2/3e ? sin 2 % , (36) 

we finally get the running parameters Mq(Q), Mjj(Q), and A t (Q). A t (Q) will be needed 
for calculating the Higgs-stop-stop and Higgs-stop-sbottom couplings. 

Bottom— sbottom sector: Here the situation is more complicated because the input 
parameter A^ is the running value at Q — m^ and all other parameters are on-shell. 
Therefore, we have to perform an iteration procedure to obtain ^((^mssm and on-shell 
m ii 2 and 9~ b . 

First we calculate the starting values for this iteration, which we denote by a super- 
script (0). From the stop sector we already know the running parameter Mq(Q). From 
Eqs. (Q) and (f|) we obtain m b (Q) SM . Together with the on-shell M^,, tan/5, and /i we 
calculate m~ and 9-* by solving the corresponding mass eigenvalue problem, see Eqs. 



1) to (5) of [ 10 1 . Next we evaluate the gluino correction term Sm^' ^ using m^, m^sM 



m~^ , and 9^ and then = m^sM + Sm^ . As the value of the running Mp is similar 

to its on-shell value we set = M5. The values tan j3 and \x are not affected by QCD 
and remain constant. The running value Mq will not be iterated because it is already 
calculated precisely enough from the stop sector. 
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The iteration procedure is as follows: 



» 



I. Tii~' J and 0g are calculated from the parameters Mq(Q), A b (Q), L> , and 



,(n-l) 



(n-l) 



2. We calculate 5mjf' n ^ according to Eq. (|28|) using ^^ 2 , 0g , and m[ n ^ (instead of 



M b ). 

» 



3. ml' = m bSM + 5m\ 

4. Sm^ 1 ] 2 and 50- are calculated from , 9~^\ and using Eqs. fl32|) and (p9|). 



5. The on-shell values m 



(») 



m-[ n )2 - 5m-r ;z , and 0f J = - 

01,2 012 oos b 



An) 2 



i(n) 



An) 



(n) 



6. SM 



(n)2 



Sm^ 2 sin 2 O-^ 1 +5m\ n> L cos 2 Q y A l) + (m 



D 6j 

(ft) _ 



( n ) 2 /)(ft) 



with 5m[ n ' ) = m[ n ' ) — M b . 



In 



bos 



-m 



2 )W sin2^ (n) ^ (n) -2M^m£ n) 



5 2 /os 



feos 6 



7. = sjMl+5Mf 2 . 



Note that all quantities without the subscript "OS" are DR running ones. 



1 



(0) 



. The iteration 
1 smaller than a 



We define an accuracy parameter A(x) by A(x) — x x (»-i)(Q) 
starts with n = 1 and stops, when A(M i = ) ), A(m fe ), and A(^) are a, 
given accuracy In our analysis we require A(x) < 1CT 5 . Thus we obtain m bi2QS , 9 bos , 
and m b (Q), which we need for the calculation of the Higgs boson decay widths. The 
consistency of this procedure is checked by calculating the on-shell from 

M l = K os sin2 6 ios + ml os cos 2 6- bos - Ml - D R (b) , (37) 
which must be equal to the input within the same order of the above accuracy. 



Decays into quarks: We first consider the tree-level as the zeroth approximation. 
As described in section [3.2| , for the decays of A and H^ 1 we take the running mass 
mt,b(jn<f>) = "^,&(?A/>)mssm in the Yukawa couplings, i.e. 

yi = -iV2a b = h b sin(3 = V2 mb< ^ mA ^ tan/3, 

v 

y 2 = -iyfta* = h t cosp = V2 m ^ m ^ cot (3 , (38) 



12 



with v = 2mw/g and the H + tb couplings, see Eq. (20) of |TD|. As outlined in 
section |3~2|, for the couplings of h° and H° to top and bottom quarks we use 



4 



'A 



m t (mh) cos a 



f sin /3 
nT>t(mH) sin a 



f sin /3 

sin (a — p) = h cos(a — p) tanp ■ 



, m m t (m H ) S M . . , m , fl m b {m H ) , . 

cos(a — p) = h sm(a — p) tanp . (39) 



One has to bear in mind that in the above procedure we have already absorbed the 
counterterms to m q with gluon and gluino exchange in the running quark mass m q (Q) 
in a*, a b , s^, and s l H . Therefore the corresponding counterterms Sa 1 (°\ 5a b (°\ Ss^ , 
and are zero[|. Only for the decays of h° and H° into bottom quarks the gluino 

contribution is just partly included. Therefore the counterterms are 

6(0) . / f3\ dTrl b 



Ss h = sin(a — (3) 



5s ^ = - cos(a - (3) . (40) 

v 

In the Higgs-squark-squark couplings G\^ k in the loop (for notation, see fllPp m q , 9 q , and 
A 9 are all taken running at the scale m^. 

Decays into squarks: The tree-level matrix element is directly proportional to the 
Higgs-squark-squark couplings G q ^ k which are functions of m q , A q , and 9 q . According 
to section f| we use running m q and A q , but on-shell 9 q . As we have absorbed the 
counterterms for m q and A q in the tree-level couplings, we have the shifts 5m q = 5A q = 
in the calculation of the one-loop corrections. The remaining counterterms are (i' ^ i 
and f ^ j) 

sol}? = - U-iyG" + (-iMJ S9 q (4i) 



for H k = {h°, H }. For H+ -> Ub* we have 

5G\% = -(-lYG i/]4 69; - (-l) j G if4 59~ b . (42) 

The tree-level coupling A q~lq 2 is independent of 9 q . Therefore no counterterms are left 
in this case. For the Higgs couplings to quarks in the loop we take Eqs. (|38D and (j39|). 



In the numerical analysis we take for the Standard Model parameters M t = 175 
GeV, M b = 5 GeV, a s (m z ) = 0.12, m z = 91.2 GeV, sin 2 9 W = 0.23, and a em (m z ) = 
1/129. 

1 Here the superscript (0) denotes the counterterms to the couplings; notation as in [[To|| . 
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Although we treat only the SUSY QCD corrections, the inclusion of the very large 
Yukawa correction to the Higgs boson sector |26| is indispensable for a realistic study. We 
thus use the formulae of ref. |27j for the radiative corrections to (m^o, m#o, a). For m#+ 
the tree-level formula is used, which is a very good approximation for m A ^> mz- 

We choose the Higgs and SUSY parameters such that they satisfy the mass bounds 
from direct searches |28| , the constraint from electroweak 5p bounds on t and b p9| using 
the formula of JJIJ, and the approximate necessary condition for the tree-level vacuum 
stability [0, 



A: 



with 



Al(Q) < 3 (M?(Q) + M? (Q) + m\) (43) 



; (Q)<3(Mj(Q)+M?(Q)+m^), 

< 2 A + m|) cos 2 /3 — \ m%, m 2 Hi = {m\ + m|) sin 2 /3 — | m|, and Q ~ M c 



m% 2 = 



(m, 



5.2 Numerical results and discussion 



In the following we show the numerical improvement of the SUSY QCD corrections to the 
widths of Higgs boson decays to quarks and squarks. First we show in Fig. [IJ the SUSY 
QCD running bottom quark mass at the scale of h°, m^irrih) (dashed line), and at the scale 
of A , m^mA) (full line), as a function of tan/3. The following values for the parameters 
are used: (M 6 (t), M , M b ) = (300,270,330) GeV, A t = 150 GeV, A b (Q) = -700 GeV 



(where Q = rrih or ttla) 



350 GeV, n = 260 GeV, and m A = 800 GeV. 




5 10 15 20 25 30 35 40 

tan (3 



Figure 1: Running bottom mass as a function of tan (3. The on-shell value is M& = 5 GeV. 
Furthermore, m(,(m/,)sM — 3.0 GeV with = 93 to 103 GeV, and m^m^sM = 2.6 GeV. 
The other parameters are given in the text. 

In Fig. |] the tree-level and corrected widths of Higgs boson decays to bottom quarks 
are shown in two ways of perturbative expansion, (i) the strict on-shell scheme (dash- 
dot-dotted line: tree-level, dash-dotted line: one-loop corrected), and (ii) the improved 
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scheme as discussed in section |372| (dashed: tree-level, full line: one-loop corrected), as a 
function of tan f3. The parameters used are the same as in Fig. |] with Q at the mass of 
the parent Higgs boson. (Notice that for case (i) A b = —700 GeV is the on-shell value.) 
In addition, the tree-level decay width using (iii) the non-supersymmetric QCD running 
mass m^m^sM is shown (dotted line). One can clearly see that the difference between 
tree-level and corrected widths decreases dramatically from (i) to (ii), especially for de- 
cays of the heavier Higgs bosons (H°, A , H^). (The curves for T(A° — > bb) are practically 
the same as those for T(H° — > bb).) In particular, the physically meaningless "negative 
width" is absent for (ii). This demonstrates the improvement of the convergence of the 
SUSY QCD perturbation expansion by the method proposed in this paper. Furthermore, 
Fig. |^a shows that (iii) is indeed a good approximation for h° — * bb in the large ijia limit, 
Eq. fl2"2"|). This is, however, not the case for the decays of the heavier Higgs particles to 
quarks, see the Figs. ||b and 0c. 

We next show the numerical improvement of the SUSY QCD corrections to the decay 
widths to squarks. Figure ^ shows the tree-level and corrected decay widths of H° — > b\b\ 
for fi = 260 GeV and -260 GeV, and H+ -> tj)\ for /j = 260 GeV, as functions of 
tan/?. The other parameters are the same as in Fig. ^. We compare the strict on-shell 
perturbation expansion of jlO) (i) and the improved one (ii) as described in section f|. 
The improvement of the perturbation expansion is again very clear. The results for 
T(A° -> &1&2) are similar to those for T(H° — > b\b\). 

We now calculate the parameter dependence of the SUSY QCD corrected widths 
T(h° -> 66), T(H° -> 66), T(H° -> bj>l), and T(A° -> bj>* 2 ) for tan/? = 30, which could 
not be presented in || [10] because the widths would have become negative. The results 
are shown in Figs. f| and § as contour plots in the (A b (Q), /x) plane for tan/3 = 30, with 
the SU(2) gaugino mass M = 117 GeV, and the other parameters being the same as in 
Fig. 0. 

As expected from Eq. (0), in Fig. ^ the corrections to T(h° — > bb) and T(H° — > bb) 
are mainly determined by fi. Note also that the \x dependence is much stronger for 
T(H° — > bb) in Fig. |]b than for F(h° — > bb) in Fig. ^|a, in accordance with the discussion 
in section 3.1. 

As can be seen in Fig. |5]a, the decay width of H° — > b\b\ becomes large for large 
|Af,(Q)| and/or large negative fi. One can see that the tree-level invariance of the decay 
width under the sign change (Ab, fi) — > (—Ab, —fi) is badly violated. This is due to the 
gluino loops. 

Finally we comment on the rrig dependence. It has been pointed out in 0, |9| that the 
gluino correction to quark modes decouples very slowly for increasing irtg (m g > rriq) with 
fixed squark parameters. This property can be understood from Eq. (|T3"|). The correction 
Ab scales as (In nig) /nig in this limit, decreasing much more slowly than the usual 1/M 2 
behavior. 
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Figure 2: Widths of Higgs particle decays into quarks as a function of tan/3. Case (i): 
dash-dot-dotted line corresponds to the on-shell tree-level, and dash-dotted to the on- 
shell one-loop result. Case (ii): dashed line corresponds to the improved tree-level, and 
full line to the improved one-loop result. Case (iii): dotted line corresponds to the tree- 
level result improved only by using rrib(Q) sm- For details, see the related text. 
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Figure 3: Widths of Higgs particle decays into squarks as a function of tan/3. Case (i): 
dash-dot-dotted line corresponds to the on-shell tree-level, and dash-dotted to the on- 
shell one-loop result. Case (ii): dashed line corresponds to the improved tree-level, and 
full line to the improved one-loop result. For details, see the related text. 
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Figure 4: Contour lines for the improved one-loop decay widths of Higgs particle decays 
into bottom quarks as a function of Ab running and \x. In (b) the light gray area shows 
the region where the on-shell one-loop result is negative. The other parameters are given 
in the text. The dark gray area is excluded by LEP and Tevatron data |E8| . 
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Figure 5: Contour lines for the improved one-loop decay widths of Higgs particles decay 
into bottom squarks as a function of Ab running and /i. The light gray areas show the 
region where the on-shell one-loop result is negative. The other parameters are given in 



the text. The dark gray area is excluded by LEP and Tevatron data p8 . 
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6 Conclusion 



The SUSY QCD corrections to the decays of MSSM Higgs bosons to quarks and squarks 
of the third gereration are phenomeno logically very important. The existing formulae for 
the one-loop SUSY QCD corrections to their decay widths have been given in the on- 
shell expansion for quarks and squarks. However, for large tan/3, these formulae suffer 
from very bad convergence of the perturbation series, which makes the numerical results 
unreliable. 

We have worked out a method for improving the SUSY QCD corrections to these 
widths. The essential point of the improvement is to define appropriate tree-level cou- 
plings of the Higgs bosons to quarks and squarks, in terms of the running quark masses 
m q (Q) both in non-SUSY and SUSY QCD, the running Higgs-squark trilinear couplings 
A g (Q), and the on-shell left-right mixing angle 8q of squarks. We have also shown numer- 
ically that this method greatly improves the corrections to the decays to bottom quarks 
and squarks. 

We finally note that the method of improvement for the Higgs boson couplings to 
quarks and squarks as presented in this paper is also useful in studying radiative cor- 
rections to other processes, such as the decays of a squark into chargino/neutralino and 



quark 



those into lighter squark and Higgs boson [IT, [T2"] , low energy four-fermi 



interactions with virtual charged Higgs boson 
bosons with quarks p|] and squarks p5 . 



33], and associated production of Higgs 
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